Search for: All records

We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicityσ_cand compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, withS₂(b) − S₂(u) ∝ l_⊥or equivalently $-E_r\propto k_{\perp }^{-2}$. HereS₂(b) andS₂(u) are second-order structure functions of magnetic field and velocity respectively. In most runs,S₂(b) develops a scaling relation $S_2\left(\mathit{b}\right)\propto l_{\perp }^{1/2}$( $E_b\propto k_{\perp }^{-3/2}$). In contrast,S₂(u) is consistently shallower thanS₂(b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initialσ_cand are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales. 

Numerical simulations have been an increasingly important tool in space physics. Here, we introduce an open-source three-dimensional compressible Hall-Magnetohydrodynamic (MHD) simulation codeLAPS(UCLA-Pseudo-Spectral,https://github.com/chenshihelio/LAPS). The code adopts a pseudo-spectral method based on Fourier Transform to evaluate spatial derivatives, and third-order explicit Runge-Kutta method for time advancement. It is parallelized using Message-Passing-Interface (MPI) with a “pencil” parallelization strategy and has very high scalability. The Expanding-Box-Model is implemented to incorporate spherical expansion effects of the solar wind. We carry out test simulations based on four classic (Hall)-MHD processes, namely, 1) incompressible Hall-MHD waves, 2) incompressible tearing mode instability, 3) Orszag-Tang vortex, and 4) parametric decay instability. The test results agree perfectly with theory predictions and results of previous studies. Given all its features,LAPSis a powerful tool for large-scale simulations of solar wind turbulence as well as other MHD and Hall-MHD processes happening in space. 

Abstract The heliosphere is permeated with highly structured solar wind originating from the Sun. One of the primary science objectives of Parker Solar Probe (PSP) is to determine the structures and dynamics of the plasma and magnetic fields at the sources of the solar wind. However, establishing the connection between in situ measurements and structures and dynamics in the solar atmosphere is challenging: most of the magnetic footpoint mapping techniques have significant uncertainties in the source localization of a plasma parcel observed in situ, and the PSP plasma measurements suffer from a limited field of view. Therefore, it lacks a universal tool to self-contextualize the in situ measurements. Here we develop a novel time series visualization method named Gaussianity Scalogram. Utilizing this method, by analyzing the magnetic magnitude data from both PSP and Ulysses, we successfully identify in situ structures that are possible remnants of solar atmospheric and magnetic structures spanning more than 7 orders of magnitude, from years to seconds, including polar and midlatitude coronal holes, as well as structures compatible with supergranulation, “jetlets” and “picoflares.” Furthermore, computer simulations of Alfvénic turbulence successfully reproduce the Gaussianization of magnetic magnitude, supporting the observed distribution. Building upon these discoveries, the Gaussianity Scalogram can help future studies to reveal the fractal-like fine structures in the solar wind time series from both PSP and a decades-old data archive. 

Abstract Parker Solar Probe observations reveal that the near-Sun space is almost filled with magnetic switchbacks (“switchbacks” hereinafter), which may be a major contributor to the heating and acceleration of solar wind. Here, for the first time, we develop an analytic model of an axisymmetric switchback with uniform magnetic field strength. In this model, three parameters control the geometry of the switchback: height (length along the background magnetic field), width (thickness along radial direction perpendicular to the background field), and the radial distance from the center of switchback to the central axis, which is a proxy of the size of the switchback along the third dimension. We carry out 3D magnetohydrodynamic simulations to investigate the dynamic evolution of the switchback. Comparing simulations conducted with compressible and incompressible codes, we verify that compressibility, i.e., parametric decay instability, is necessary for destabilizing the switchback. Our simulations also reveal that the geometry of the switchback significantly affects how fast the switchback destabilizes. The most stable switchbacks are 2D-like (planar) structures with large aspect ratios (length to width), consistent with the observations. We show that when plasma beta (β) is smaller than one, the switchback is more stable asβincreases. However, whenβis greater than 1, the switchback becomes very unstable as the pattern of the growing compressive fluctuations changes. Our results may explain some of the observational features of switchbacks, including the large aspect ratios and nearly constant occurrence rates in the inner heliosphere. 

Abstract A variety of magnetosphere‐ionosphere current systems and waves have been linked to geomagnetic disturbance (GMD) and geomagnetically induced currents (GIC). However, since many location‐specific factors control GMD and GIC intensity, it is often unclear what mechanisms generate the largest GMD and GIC in different locations. We address this challenge through analysis of multi‐satellite measurements and globally distributed magnetometer and GIC measurements. We find embedded within the magnetic cloud of the 23–24 April 2023 coronal mass ejection (CME) storm there was a global scale density pulse lasting for 10–20 min with compression ratio of . It caused substantial dayside displacements of the bow shock and magnetopause, changes of and , respectively, which in turn caused large amplitude GMD in the magnetosphere and on the ground across a wide local time range. At the time this global GMD was observed, GIC measured in New Zealand, Finland, Canada, and the United States were observed. The GIC were comparable (within factors of 2–2.5) to the largest ever recorded during 14 year monitoring intervals in New Zealand and Finland and represented 2‐year maxima in the United States during a period with several Kp7 geomagnetic storms. Additionally, the GIC measurements in the USA and other mid‐latitude locations exhibited wave‐like fluctuations with 1–2 min period. This work suggests that large density pulses in CME should be considered an important driver of large amplitude, global GMD and among the largest GIC at mid‐latitude locations, and that sampling intervals are required to capture these GMD/GIC. 

Abstract This Letter reports the first observation of the onset of fully developed turbulence in the solar corona. Long time series of white-light coronal images, acquired by Metis aboard Solar Orbiter at 2 minutes cadence and spanning about 10 hr, were studied to gain insight into the statistical properties of fluctuations in the density of the coronal plasma in the time domain. From pixel-by-pixel spectral frequency analysis in the whole Metis field of view, the scaling exponents of plasma fluctuations were derived. The results show that, over timescales ranging from 1 to 10 hr and corresponding to the photospheric mesogranulation-driven dynamics, the density spectra become shallower moving away from the Sun, resembling a Kolmogorov-like spectrum at 3R_⊙. According to the latest observation and interpretive work, the observed 5/3 scaling law for density fluctuations is indicative of the onset of fully developed turbulence in the corona. Metis observation-based evidence for a Kolmogorov turbulent form of the fluctuating density spectrum casts light on the evolution of 2D turbulence in the early stages of its upward transport from the low corona. 

Abstract This paper addresses the first direct investigation of the energy budget in the solar corona. Exploiting joint observations of the same coronal plasma by Parker Solar Probe and the Metis coronagraph aboard Solar Orbiter and the conserved equations for mass, magnetic flux, and wave action, we estimate the values of all terms comprising the total energy flux of the proton component of the slow solar wind from 6.3 to 13.3 R ⊙ . For distances from the Sun to less than 7 R ⊙ , we find that the primary source of solar wind energy is magnetic fluctuations including Alfvén waves. As the plasma flows away from the low corona, magnetic energy is gradually converted into kinetic energy, which dominates the total energy flux at heights above 7 R ⊙ . It is found too that the electric potential energy flux plays an important role in accelerating the solar wind only at altitudes below 6 R ⊙ , while enthalpy and heat fluxes only become important at even lower heights. The results finally show that energy equipartition does not exist in the solar corona.